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The frequency spectrum of density fluctuations is calculated for a fluid whose molecules possess 
an internal degree of freedom which is weakly coupled to the translational degree of freedom of the 
fluid. Irreversible thermodynamics is used to obtain an equation of motion for the internal degree of 
freedom. This equation plus the linearized hydrodynamic equations are solved for the frequency spec- 
trum of density fluctuations. The results are compared with a similar calculation involving a frequency 
dependent volume viscosity. The results are identical for structural relaxation but there is a difference 
for thermal relaxation. The origin of the difference is discussed and the magnitude of the difference 
is examined for CCLj and lor CS2. 
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1. Introduction 

The spectrum of light scattered by density fluctua- 
tions in a fluid is proportional to frequency spectrum 
of density fluctuations [l]. 1 The development of the 
gas laser has made it possible to measure the spectra 
of long wavelength density fluctuations and thereby 
study fluids in the frequency region of a few giga- 
cycles down to essentially zero cycles/sec. 

For dense fluids this frequency region corresponds 
to slowly varying processes whose time dependence 
may be reasonably expected to be described by the 
equations of irreversible thermodynamics. When 
considering density fluctuations the appropriate 
equations are the linearized equations of hydrody- 
namics plus two additional relations among the vari- 
ables so that a solution to the initial value (or boundary 
value) problem is possible. 

In this paper we shall consider density fluctuations 
in a fluid whose molecules possess an internal degree 
of freedom which relaxes toward local equilibrium with 
the density and/or the temperature [2, 3]. This is carried 
out in section 2. In section 3 this calculation is com- 
pared with an earlier calculation [1] where the internal 
degree of freedom was assumed to result in a frequency 
dependent volume viscosity. The introduction of a 
frequency dependent volume viscosity is based on 
statistical mechanical arguments [4, 5]. 



The comparison is made in two limiting cases. In 
the first case it is assumed that the internal variable 
depends only on the local temperature (thermal 
relaxation). It is found that the two calculations yield 
similar, but not identical, results. The origin and magni- 
tude of the difference are examined. Carbon tetra- 
chloride and carbon disulfide are used to illustrate the 
discussion. In the second case it is assumed that the 
internal variable depends only on the local density 
(structural relaxation). The two calculations yield 
identical results. 

No comparison is made when both thermal relaxa- 
tion and structural relaxation are present. 



2. Calculation 

The thermodynamic statement relating 
eters for our system is [6] 



dU-- 



TdS+^dp+Azdt 



the param- 



(i) 



where U is the energy, S is the entropy, T is the 
temperature, p is the density, p<> is the equilibrium 
value of the density, p is the pressure, £ represents an 
internal degree of freedom and Ag is the partial deriva- 
tive of the Helmholtz free energy with respect to £, 



At = (dAld£) p . T . 



(2) 



*TW8 work was sponsored by the Advanced Research Projects Agency of the Department We shall further require that At = when the internal 
hi Defense. » ' 

1 Figures in brackets indicate the literature references at the end of this paper. degree of freedom is ill local equilibrium Wit ll the 
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density and the temperature (whether or not the dens- 
ity and temperature have reached equilibrium values). 
This may be expressed as 
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(13) 



Z = Z(P,T) 



(4) 



where r— (LA^) l . 

A similar equation is obtained for the entropy. 



is the local equilibrium value of £. 

A kinetic equation for £ is obtained by applying the 
methods of irreversible thermodynamics to the sys- 
tem with the result 



d^/dt- 



-LA* 



(5) 



where the kinetic coefficient L is > 0. As we are con- 
cerned with small deviations from equilibrium, A% is 
expanded to obtain 

^=.-L[At P (p- Po )+A ( AT-To)+AK(t-Z )] (6) 
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Equations (8), (9), (10), (13), and (14) constitute a 
complete set of equations. The next step is to obtain 
the Fourier-Laplace (space-time) transform of this 
where A^= (dA € ldp) T ;t, etc. If we expand eq (3) and set of equations. The result is 
substitute this in eq (6) we obtain 

zpk(z) + potykU) — p k , (15) 

d€ldt = -LA H (€-€). (7) 



The equations of motion for the system are the lin- (z-hbok 2 )\bk(z) 
earized hydrodynamic equations: 
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Here i// = div v, r} 8 is the shear viscosity, rj v is the vol- an( * 
ume viscosity and k is the thermal conductivity and 
m is the mass of a molecule. We relate the pressure 
and entropy to the density, temperature and f by the 
thermodynamic relations 
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Here we have made use of thermodynamic relation- 
The internal variable £ is eliminated by combining ships which are developed in the appendix. The trans- 



eqs (6), (11), and (12). For the pressure we find 
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formed quantities are 

friz) = J*dte-* j A«*-*[p(r, t) -p,]. (18) 
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7a- (z) and i|/a-(z) are obtained in the same way. In eq 
/4 



(16) bu^ ( -j) s + r) r \ Impu, Co is the low-frequency 
sound speed, y is the ratio of specific heats. 



G(z)=Z 4 T(l-C l IC r ) 

+ z 3 [l + ak 2 T+b»k 2 T(l-C l IC r )] 

+ z 2 [ak 2 + b k 2 + ab k 4 r + C 2 Fr ( 1 - C//C„ ) ] 

+ z[abok 4 + Clk 2 (l + ak 2 Tly)] 

+ ak 2 Clk 2 ly. (23) 



/3r is the isothermal coefficient of thermal expansion, 
tp=-AeplAff and ijf=-Af T IAtf. In eq (17) C,=T {) A^ 2 
and ^/ = K/poCr. The quantity C/ represents the con- 
tribution of the internal degree of freedom to the hor purposes of comparison let us recall the corre- 
specific heat when the density is held constant. sponding expressions derived earlier for the case of a 

The right-hand sides of eqs (16) and (17) may be frequency dependent volume viscosity [1]. There 
simplified by noting that in equilibrium 



-^-dp=^-dp + ^^dT 

mp„ yp, y 



and 
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This simplification is permissible because we are 
interested in computing the correlation function 
(pkP-k(z)) in terms of the equilibrium correlation 
function ( |pa-| 2 ). 

To obtain (pkp a-(z)) we solve eqs (15-17) forpA-(z) 
in terms of the initial values p*., 7/- and i//a-. Only the 
terms involving pi- are of interest because 

<PA.r_A> = <PA^A) = 0. 

3. Comparison With a Calculation Using 
a Frequency Dependent, Single Relax- 
ation Time, Volume Viscosity 

3.1. Thermal Relaxation 

We are interested in the case where the relaxing 
parameter is independent of the density, i.e., ij p = 0. expressions are not equivalent. This may be realized 
The correlation function is found to be 



F(z)=z :i; r + z 2 [l+T(ak 2 + b {) k 2 )] 

+ z[ak 2 + boh? + b x k 2 + abok 4 r + Clk 2 r( 1 - 1/y) ] 
+ abok 4 + abik 4 + Qk 2 ( 1 - 1 /y ) . (24) 

and 

G(z) = rz A + z :i [ 1 + r ( ak 2 + bole 1 ) J 

+ z 2 [ak 2 + bok 2 + b t k 2 + r(Clk 2 + abok 4 ) ] 

+ z[Qk 2 + aW 4 + ab x k 4 + fCj# 2 aA: 2 /y] 

+ C 2 ,A 2 a£ 2 /y. (25) 

[Quantities with a tilde refer to the calculation in ref. 1 . | 
Here 6<> is the nonrelaxing part of the kinematic vis- 
cosity, the relaxing part is b\l(\ + Iojt) 
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where 
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Direct comparison of 6' and G shows that the two 
tpressions are not equivalent. This may 
by noting that G contains the combination 

z^bok 2 T{\-CilC v ) +£ 1 C 2 k 2 {\-C l IC l> )r 
while the corresponding terms in G are 

Z*bok 2 T + Z 2 Qk 2 T. 

If we consider the limiting case where the internal 
degree of freedom provides the only loss mechanism 
for translational degrees of freedom it is possible to 
make G and G into equivalent expressions. In this limit 
where a = b { ) = 

G(z)=z 4 T(l-C,IC r )+z :i 

+ *Qh*T(\-CtlC p ) +zCW (27) 



and 
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G(z)=-z 4 f+z 3 
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A simple comparison shows that if 

t = t(\-C,IC v ) 



— ) + zClk 2 . (28) 
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(34) 



(29) 



then 6' = 6' in this limit. It also can be shown that F = F 
under these conditions. 

In an experiment what is measured is not F(z)/G(z) 
but the real part of FjG when z— ico: 



(pkp-kM) _ n F(ico) 

7T T\ — Wx — Ke • 

G(ico) 



<l/*l 



(30) 



D 2 = -co :i [l+ak 2 T + b {) k 2 T] 

) C 2 k 2 (\ + ak 2L ) + ab {) k 4 



+ ak 2 TClk lC ^(\-- 

C r \ y, 



Examination of these expressions reveals; (a) they are 
not identical term for term, (b) there will be no sig- 

A comparison of the complete expressions requires nificant difference in their values if 

that we compare 



G ( iio ) 
To do this we write 
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(32) 



G(m) D 2 + D 2 
From eqs (22) and (23) it follows that 

/V, - - co 2 [ 1 + ak 2 r + M- 2 r] + «/U 4 -f Cg£ 2 (l -1/y) 
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Here use has been made of eq (29). 

The corresponding expressions for N h etc., are 

yV, =- co 2 [l + ak 2 T+ b {) k 2 r] + ai»A; 4 

+ Co 2 A 2 (l-l/y)+C ( U' 2 aA- 2 r^(l-l/y) 



a = ak 2 rCilC r 



(35) 



is small compared to 1. 

Using the parameters listed in ref. 1. we find that a 
is 6.4 X10- 4 for CC1 4 at 20 °C and is 2.0 X 10~ 2 for 
CS 2 at 20 °C. The wavevector was taken to be 10 5 cm -1 . 
This could be increased to a = 0.2 for CS2 if we con- 
sider scattering angles on the order of 180° using the 
He — Ne laser as the light source. 

The spectra predicted by eqs (31) and (32) have 
been evaluated for CCI4 and CS2 using the parameters 
listed in ref. 1. For CCI4 there is no difference between 
the two curves. This is not surprising considering the 
small value of a. 

The two curves are not identical for 

CS 2 (A: = 3Xl0 5 cm- 1 ). 
The central components are shown in figure 1. The 




2 3 

oj.IO 8 rod. /sec 



FIGURE 1. The central component of light scattered from CS 2 as 
predicted by eq (32), curve A and as predicted by eq {3D curve B. 

The parameters are taken from table of ref. 1. The wave vector /.=3 X ]()■> ( - m -'. 
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curve labeled A was obtained using a frequency de- 
pendent volume viscosity (eq 32). The curve labeled B 
was obtained using the results of this calculation 
(eq 31). The difference in the amplitude at co = is 
due entirely to the term 



c 



(Wa&T-^ (1-1/y) 



in N \. Such a quantity is missing in N\. It occurs as a 
result of multiplying the thermal diffusivity and the 
frequency dependent volume viscosity. We note for 
cd = 0, the present calculation yields the same ampli- 
tude as the theory which does not include internal 
relaxation processes. 

The Brillouin components for CS2 are shown in 
figure 2. The differences are due to the relatively 
large value of a and to the different values of ojht 
where o>/< is the Brillouin frequency. For curve A, 
co#f=78 while for curve B, o)ht= 113. The large value 
of a is the dominant feature as u)n depends most 
strongly on k when q)hT > 1. 

The possibility of finding a situation where a is 
appreciably greater than 0.2 appears to be remote. 
Longer relaxation times exist for rotational isomers 
[3] , however, Ci/C v is very small so a is still much 
less than unity. 

3.2. Structural Relaxation 

Next let us examine briefly the other limiting case 
which involves a structural relaxation rather than a 
thermal relaxation. A purely structural relaxation 
occurs when £7 = 0. The examination is speedily con- 
cluded by observing that if we formally identify A in 
eq (16) with bi/T = C 2 x — C§ (see eq (22) of ref. 1) then 
this case and the frequency dependent volume vis- 
cosity calculation yield identical results. Also no 
scaling of the relaxation times is needed. 
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4. Summary and Discussion 

The frequency spectrum of density fluctuations in a 
fluid containing an internal degree of freedom has 
been calculated. The procedure used was to obtain an 
equation of motion for the internal variable £ by the 
methods of irreversible thermodynamics and then to 
solve the resulting set of coupled equations for the 
spectrum of the density fluctuations. 

This calculation was compared with an earlier 
calculation which treated the internal degree of 
freedom as a frequency dependent volume viscosity 
rather than as a separate variable with an equation 
of motion of its own. Two limiting cases were examined. 
The structural relaxation case, fr = was found to 
yield identical results when compared with the 
earlier calculation. The thermal relaxation case, 
f P = does not agree exactly with the earlier calcula- 
tion. The difference results from the ways the stresses 
due to the internal variable are treated. In the present 
calculation the stresses are additive only to first order 
[7|. In the earlier calculation they were strictly additive. 

The equation of motion approach worked out in 
this paper has an intuitive appeal in that one can see 
directly how the internal variable enters the problem 
and how it is related to the temperature and the 
density. The introduction of a frequency dependent 
volume viscosity in the case of thermal relaxation is 
required by statistical mechanical considerations 
[4, 5] and is therefore the preferred way to handle 
internal degrees of freedom which are weakly coupled 
to the translational degrees of freedom. 

This calculation could be adapted to a chemically 
reacting system. This aspect of fluctuation theory and 
light scattering is discussed at length in recent publica- 
tions [8, 9| so we shall not discuss it here. 



5. Appendix 

Starting with eq (1) the differential expression for 
the Helmholtz free energy is 



dA^SdT+^dp + Afdg. 
Po 



(36) 



The Maxwell relations implied by eq (36) are 
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FIGURE 2. The Brillouin components of light scattered from CS2 a 
predicted by eq i32), curve A and as predicted by eq 131) curve B. 

The parameters taken from table of ref. 1. The wave vector /, =3 x 10"' cm '. The hor 
zontal lines indicate the full width at half maximum for the components. 
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In equilibrium A* = so that Finally we note that 

dAt = A {p dp+Ae T dT+AMd€ = 0. mQ = /dp 

J \dpjr 



It follows that 



mp p T Cl_ /dp 



At 9 = -Au€p \ d T/ p 

and 

AtT=-Au£T ( d A\ = Cp-Cv 

\dpjr 



where £ P = (d£/dp)r and £ T = (d£ldT) p . 

If we combine these results with eqs (11) and (12) we 

obtain 



\dp/ T poTfir 



A most interesting discussion with C. I. A. Stegeman 



dp/r.€ \dp/T 
dp\ _ (dp 



'a„\ /n.A provided the stimulus for this work. 
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